Field-induced metal-insulator transition in the c— axis resistivity of graphite 
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We show that the resistivity perpendicular pc and parallel pa to the basal planes of different 
graphite samples show similar magnetic- field-driven metal-insulator-transitions at a field _Bc ~ O.f T 
applied parallel to the c— axis. Our results demonstrate the universality of the recently found 
scaling in pa of graphite and indicate that the metallic-like temperature dependence of pc is directly 
correlated to that of pa- The similar magnetoresistance found for both resistivities, the violation 
of Kohler's rule and the field-induced transition indicate that the semiclassical transport theory is 
inadecuate to understand the transport properties of graphite. 



Although a considerable amount of studies has been 
performed on graphite and related compounds, their 
transport properties are still not well understood. The 
scientific interest on graphite has been recently renewed 
by new magnetization [0 and transport |Q results on 
highly oriented pyrolytic graphite (HOPG). These show 
irreversible magnetization |0] that, upon sample, its pre- 
vious thermal treatment and magnetic field orientation, 
resembles that of a superconducting loop even at room 
temperature suggesting the existence of localized su- 
perconducting domains at high temperatures in topo- 
logical disordered regions This result is indirectly 
supported by the recently found magnetic-field-driven 
superconductor-insulator-type transition (SIT) in the in- 
plane resistivity of a HOPG sample j^. Remarkably, this 
field-driven transition showed a similar scaling as found 
for two-dimensional (2D) disordered superconductors as 
well as for Si MOSFETs g. 

It has been speculated that in graphite impurity- 
assisted hopping mechanism may affect the electronic 
transport along the c— axis. In this case the conduction 
across the graphite planes might be influenced by the in- 
plane transport in a given temperature (T) range. The 
aim of our research is to check whether the c— axis resis- 
tivity Pc{T, B) shows a similar SIT with a similar scaling 
as the in-plane resistivity pa{T,B). Evidence for this 
speculation has important consequences to solve an old 
and basic problem of the transport in graphite, namely, 
to what extent the metallic-like behavior of pc is an in- 
trinsic property of the transport along the c— axis and 
whether the semiclassical theories are appropriate to ac- 
count for the observed effects. Furthermore, it is impor- 
tant to prove experimentally the universality of the 2D 
scaling found in graphite; for this aim different graphite 
samples with different resistivities have been measured 
in this work. 

We have studied the resistivities Pa{T, B) and Pc{T, B) 
of four samples: two HOPG samples, sample 1 from Ad- 



vanced Ceramics (rocking curve FWHM ~ 0.4°), sam- 
ple 2 (from the Research Institute "Graphite" (Moscow) , 
FWHM ~ 0.4°) and samples 3 and 4, Kish graphite crys- 
tals. The samples typical length and width were 2 mm 
and thickness between 0.1 mm and 0.5 mm. The in-plane 
resistivity pa was measured attaching four contacts on 
the sample surface with silver paint. The absolute val- 
ues were obtained with other configuration that assures a 
uniform current distribution. The c— axis resistance was 
measured attaching two contacts on each of the surfaces 
(one in the middle and other surroundig it and contacted 
on the rest surface). A resistance bridge (LR700) was 
used with an ac current < 1 mA. The field was applied 
normal to the graphene layers. The temperature stability 
was better than 2 mK in the whole range. 

The left panels of Figs. 1 to 3 show both resistivities 
as a function of T at constant fields for samples 1 to 3, 
respectively. At zero applied field the basal-plane resis- 
tivity shows the typical behavior reported in the litera- 
ture: whereas for sample 1 pa shows a metallic behavior 
up to T - 100 K, [r^ 220 K for sample 2), sample 3 - 
which shows the lowest resistivity ~ up to the highest 
measured T ^ 270 K. Similar absolute values as well as 
the observed maximum of pc for samples 1 and 3 were 
reported in the literature ||, |^ . Note also that sample 2 
with the highest pc(0) shows no maximum but a satura- 
tion below ^ 10 K. The basal-plane rest resistance ratio 
RRR = R{mOK) / R{A.2) has been tentatively taken in 
the literature as a measure of sample perfection H]. A 
speculative interpretation of the observed non-monotonic 
behavior of pc vs. RRR has been given in terms of the in- 
fluence of laminar defects and microcraks. Samples with 
RRR < 10, for example, were considered of poor per- 
fection. We note however, that this sample classification 
might be misleading because topological defects may play 
a role other than scattering centers. It has been recently 
demonstrated that defects can lead to a local enhance- 
ment of the density of states and generate superconduct- 
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ing as well as ferromagnetic domains ||3||. In this case 
large RRR values do not necessarily mean high quality. 
We believe that a broad spectrum of transport data of 
graphite may need to be reinterpreted. 

A clear magnetic field-driven transition from a 
metallic- to a semiconducting-like behavior is observed 
in both Pa and pc, see Figs. 1 to 3. Remarkably, this 
transition is observed at fields i? 0.1 T in both resis- 
tivities. The observed behavior in pa{T, B), in particular 
the leveling at low T, is similar to that reported in 2D 
films §. The minimum in Pa{T), that develops at inter- 
mediate fields in sample 3, is similar to that measured in 
a Si MOSFET §] but located at 100 times higher T and 
observed at 10 times smaller fields §. In spite of these 
similarities, the transition in graphite may have other ori- 
gin as in Si MOSFETs. We note that the metal-insulator 
transition (MIT) in 2D electron systems is measured in 
the in-plane resistivity as a function of the magnetic field 
applied parallel as well as normal to the film plane E|. 
Recent high-resolution angle-dependent experiments Q 
reveal however, that the MIT in graphite is driven mainly 
by the perpendicular component of the field providing a 
clear evidence of the 2D character of graphite in agree- 
ment with recent theoretical work ||ll|. 

In the ideal case the separatrix line obtained at a crit- 
ical field B — Be is a line where the resistivity remains 
temperature independent and separates the metallic from 
the semiconducting branch. It has been shown that for 
2D disordered superconductors as well as for Si MOS- 
FET's, for fields near the MIT the in-plane resistivity 
follows a scaling when it is plotted as function of the 
scaling variable \B — Bc\T~^^°', being a a critical expo- 
nent. In a similar way we determined Be as well as a from 
PaiT, B), see the right upper panels in Figs. 1 to 3. Tak- 
ing into account that the separatrix line is not an ideal 
horizontal line in our samples, we obtain a reasonable 
scaling in a restricted temperature range. Both parame- 
ters are similar for the three samples and agree with those 
reported in Q. The exponent a = 0.6 ± 0.05 coincides 
with that for Si MOSFETs ^ and uhrathin a-Bi films 
p3| . A similar scaling with the same exponent a but 
with a smaller Be is obtained for pc at 30 K < T < 60 K 
for samples 1 and 3. Sample 2 c— axis resistivity shows no 
metallic behavior at zero field and therefore no apparent 
scaling. Based on this scaling a SIT in MOSFETs has 
been suggested [|4|. The evidence for Cooper pairing of 
the carriers, however, is still a matter of discussion in the 
literature and no consent has been achieved yet jH . 

As shown in superconducting amorphous In-0 films 
[ p5[ the scaling can be considerably improved if one as- 
sumes that in addition to the contribution to the conduc- 
tivity that shows the SIT, an additional T— dependent 
contribution affects the resistivity. It has been speculated 
that this contribution can be due to the conductivity of 
normal electrons p5[ | . Within the same idea and assum- 
ing that the intrinsic c— axis resistivity of graphite has 



a non-metallic T— dependence which adds as background 
to the field dependent part, we show in the low right 
panels of Figs. 1 to 3 the scaling obtained if we subtract 
from the measured Pc(T, B) the line Pc{T, Be) or a linear 
T— dependence ap^T (sample 2, a is a free parameter and 
Po = /Oc(r ^ 0, S - 0.1 T)). In this case the T-range for 
scaling is considerable enhanced. Obviously, a similar in- 
provement of the T— range for scaling is obtained also for 
the in-plane resistivity after subtracting a T— dependent 
background. The value of the exponent a remains the 
same. These results indicate the universality of the MIT 
found in graphite for both resistivities and of its scaling. 

Graphite can be considered as a dilute electron system 
with a Coulomb coupling constant = l/(7rn2_D)^/^a|j 
where the 2D carrier density is given by n2D — nsod with 
d = 0.335 nm the interplane distance, = eh^ /e^m* the 
effective Bohr radius, e = 2.8 the dielectric constant and 
TO* the effective mass of the carriers. Taking literature 
values for the majority (minority) carriers | [l6t one ob- 
tains rs ~ 5(10) suggesting that the 2D electron system 
in graphite might be considered as a strongly correlated 
liquid and the effects due to the electron-electron interac- 
tion should be important. We note that n2D ~ 10^^ cm~^ 
practically coincides with that reported for a Si MOSFET 

i- 

One important characteristic observed in several of 
the measured 2D systems with a MIT is that the crit- 
ical sheet resistance Rg at the separatrix is of the order 
of the quantum unit h/e"^. Within the scaling theory 
of a SIT, R, ^ h/Ae^ ~ 6.4 kfi at B^. Assum- 
ing that the input current decreases exponentially in the 
c— direction with a penetration depth A and taking into 
account the distance d between graphene sheets, we esti- 
mate Rs — {b/l)R/{l — exp(— d/A)) for a single graphene 
layer, where R is the measured longitudinal resistance 
at the sample surface, b the sample width and / the dis- 
tance between voltage contacts iQ . With measured pa- 
rameters we obtain at the separatrix Rg ~ a(/i/4e^) with 
a = 1/3.2(1/3.7) for sample 1 (2), whereas for sample 3 
a — 1/27, using A 10 fim obtained from measurements 
with different configurations. The difference in sheet re- 
sistance between samples 1,2 and 3 may be related to a 
difference in the coupling between layers due to lattice 
defects, like stacking faults or in-plane defects, or may 
indicate the lack of an unique critical resistance as ob- 
served in disordered thin films ||l^ . 

We note that lattice defects may not only affect the 
transport as scattering centers but they may contribute 
to enhance the electronic coupling between layers giv- 
ing rise to a quasi 3D electronic spectrum with coher- 
ent transport along the c— axis. Field-angle dependent 
experiments around the direction parallel to the planes 
show a weak coherent peak in the c— axis resistivity for 
Kish graphite (but not for the HOPG samples) as ex- 
pected from theoretical considerations for coherent trans- 
port Ip^]. Regarding defects within the graphene layers 
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that may generate superconducting domains |^ , we note 
that for sample 3 the metaUic behavior of pc measured 
at low T remains at all fields irrespective of the behavior 
of Pa which shows an insulating-like behavior at fields 
> 300 Oe. This behavior in Kish graphite might be due 
to a "phase separation" , i.e. to the presence of both in- 
sulating and metallic domains within the planes. 

The T— dependence of pc in HOPG has been inter- 
preted in terms of a semiclassical model with pa- 
rameters like the carrier density n(T), the stacking-fault 
spacing / and the potential barrier formed on the plane 
of stacking disorder AE. The increase of pc with T at 
T < 100 K was attributed to an interlayer transfer of car- 
riers which provides the mixing of the c-axis and in-plane 
transport. It was further speculated that the interaction 
of the carriers with phonons on the planes determines 
the T— dependence d). We note first that according to 
the model ^ and from the comparison of the absolute 
values of pc of our samples, sample 2 should have an ex- 
tremely large stacking-fault density. This is at odd with 
x-ray diffraction results that indicate that this kind of 
lattice defects is absent in high quality HOPG (2l]. Sec- 
ond, from our results it is obvious that the increase of 
Pc with r at T < 100 K cannot be attributed simply to 
electron-phonon scattering. Nevertheless, the concept of 
conduction-path mixing mechanism is appealing. This 
can be realized if the conduction-path of the electrons 
along the c— axis is in part short circuited by lattice de- 
fects or impurities as was usually speculated in litera- 
ture. However, the apparent similarity of the behavior 
obtained in the three samples may indicate an intrinsic 
origin for the conduction-path mixing. A situation simi- 
lar to that found in layered high-temperature supercon- 
ductors, where in-plane charge fluctuations influence the 
c— axis transport may be possible in graphite. Our 
results cast also doubts about the estimate of ~ 0.3 eV 
for the interlayer transfer integral used in literature [p^ . 
This is not surprising since neither the electron-electron 
interaction nor charge fluctuations (23| were taken into 
account. 

Regarding the field dependence of pc one may argue 
that the field can decrease the carrier concentration at 
low T inducing an increase in pa and pc- However, this 
appears to be at odd with the increase of n{T) by a fac- 
tor of two with a field of 0.2 T obtained from magnetore- 
sistance measurements and within a semiclasical model 
p^ . As noted in the literature |Q the mean mobil- 
ity in graphite cannot be estimated rigurously with the 
usual two-band model since the resistivity does not fol- 
low the predicted dependence. Arbitrary definitions 
to circumvent this problem did not provide a satisfac- 
tory description of the measurements. Our results indi- 
cate clearly that Kohler's rule is not fulfilled in graphite, 
neither in the a nor in the c— direction -as can be seen 
in Fig. 4(a) for pa- implying that the semiclassical pic- 
ture does not apply. Furthermore, we demonstrate in 



Fig. 4(b) that the low-temperature magnetoresistance in 
the c-direction, shown for two Kish graphite samples and 
sample 1, is basically the same as the in-plane one. This 
result also indicates the correlation between both con- 
ductivities and contradicts the semiclassical model since 
for i?| |c| I J no magnetoresistance in the c— direction is ex- 
pected. 

Concluding, we proved that both in-plane and out-of- 
plane resistivities in graphite show a MIT with a scal- 
ing similar to 2D electron systems and indicate that the 
metallic-like behavior of pc{T,B) is related to pa{T,B). 
Our overall results show that semiclassical models cannot 
account for the magnetotransport properties of graphite. 
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FIG. 1: Left panels: In-plane pa and c— axis pc resistivities as a function of temperature at constant fields for sample 1. The field 
step between two curves is 0.01 T. Right panels: Scaling behavior of the resistivities. The experimental line pc{Bc = 0.04 T) 
has been subtracted from the pc-data. The symbols are the data points obtained between 30 K and 120 K. The symbols in the 
lower figure correspond to the data taken between 40 K and 270 K. 
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FIG. 2: The same as Fig. 1 but for sample 2. Right pannel: The symbols are the data points obtained between 30 K and 
200 K. The symbols of the lower figure are obtained from the data between 30 K and 110 K. The parameter a = —0.0013 K"'^ 
and po = 4.7 mfim. 




FIG. 3: The same as Fig. 1 for sample 3, pa (2K,0) = O.OOlS/Ltfim. Right pannel: The symbols are the data points obtained 
between 50 K and 130 K. The symbols of the lower figure correspond to the data taken between 90 K and 200 K. 
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FIG. 4: (a) Kohler's plot of the in-plane magnetoresistance data for sample 1. Between 10 K and 270 K each curve has been 
taken every 10 K. (b) Magnetoresistance Ap = pa,c{B) — pa.c{0) / Pa.ciO) as a function of applied field at T = 2 K for samples 
1, 3 and 4. 



